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, Primary Cosmic Ray Chemical Composition is investigated in energy re- 

gion close to 10^^ eV. Studies are based on comparisons of 7-hadron fam- 

o 

a^ 
o , 

\^ ilies generated by means of quasi-scaling model MCO. It is shown, that 

o 

all characteristics of observed families, including their intensity, are in a 
^ ' very good agreement with simulated event properties at the normal chemi- 

V ' 

Qh' cal composition and are in disagreement at heavy dominant compositions. 



ilies observed by Pamir and Pamir-Chacaltaya Collaboration, with fam- 



Code CORSICA with VENUS and DPM models also contradicts with 
experimental data of families. One- and multi-dimensional methods of 



X 

I recognition of Fe-like families is worked up and approved. They are based 

on family characteristics sensitive to atomic number of induced nuclei and 
are not correlated between each others. It is shown that the fraction of 
Fe-like families is consistent with the normal chemical composition and 
strongly contradicts to heavy dominant ones. The success of MCO model, 
in description of families properties, is due to large inelasticity coefficient 
of soft interactions at superhigh energies. 
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1. Introduction 

The investigation of Primary Cosmic Ray (PGR) Chemical Composition(CC) is one of the 
key problem for understanding of PGR origin, properties of radiation sources and interstellar 
or extra Galaxy space, which PGR passes from a source to the Earth. GG is well investigated 
up to energy lO^^eV. By direct methods GG is studied up to energy 10^^ eV H . However 
the reliability of results of these work is not great, since they base on small statistics. At 
larger energies sources of the information about GG is Extensive Air Showers (EAS) or 
families of 7-quanta and hadrons, registered by X-ray emulsion chambers . Despite of long- 
term researches the results are very inconsistent. The statements concerning GG swing from 
so-called normal |^ and even with proton dominant to heavy ||^, |^ and superheavy |^ 
compositions. The appropriate data with the references are brought in Table |. 

Whereas the knowledge of GG in the region 10^^ - 10^^ eV is especially important as a 
complete spectrum of PGR in this area has a "knee". 

The reason could be either spectrum bend of PGR (due to decrement some the component 
from PGR content) or change properties of inelastic interaction between cosmic ray particles 
and Earth's atmosphere nucleus at these energies [|l], 

The present work is devoted to research of PGR GG in the energy region directly after 
bend of a energy spectrum of PGR. The 7-hadron families are registered and processed by 
Pamir and Pamir-Ghacaltaya Gollaborations. The characteristic of the families are used for 
this task. 

Running ahead we shall notice, that efficiency of families generated by a nucleus, e^, 
strongly depends on their atomic number decreasing with increase of A. Therefore the chem- 
ical composition of families differs much from GG of PGR: 

/a = eA X CAin^A X Ca) (1) 

For comparison see Table | and |I[ 

Table |I| explains the main difficulty of investigation of PGR GG based on data about 
families. Even if characteristics of families induced by heavy nuclei strongly differ from that 
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TABLE I: Used chemical composition of PGR at Eq = 10^^ eV , Ca 



Composition 


P 


He 


CNO 


SiMg 


Fe 


{In A) 


Normal |, || 


40 


20 


10 


10 


20 


1.7 


Heavy |, § 


15 


10 


17 





58 


3.0 


Superheavy 


7 


5 


12 


6 


70 


3.4 



TABLE IL Chemical composition of families, /^% 



Composition 


P 


He 


CNO 


SiMg 


Fe 


Normal [|, | 


75. 


16.6 


3.2 


2.4 


2.8 


Heavy |, § 


56. 


16.6 


11. 





16.4 


SuperveavyQ 


40. 


12.6 


12. 


4.4 


31. 



of proton families, influence of the former to the average characteristics of families is not 
great, since CC of families is strongly enriched by protons. Even at heavy composition of 
PGR, in which iron makes about 60%, they provide only 16.4% of families. How Table |I| is 
received and what is possible to make for the analysis of PGR GG by data on families will 
be discussed in the subsequent sections. 

A 7-hadron family is a result of Nuclear Electromagnetic Gascade (NEG) developed in 
the Atmosphere after interaction of PGR particles with nuclei of air in the top part of 
the Atmosphere. This process is very complex, multistep and branching. Therefore to 
get information from characteristics of families one needs a comparison with Monte-Garlo 
calculation based on some mode including both nuclear and electromagnetic cascades. 

Not less than ten of such models [^, |10|, |Tl], |T2|, [13|, |T5|, [16|, were developed up 



to day. All recent models of nuclear interactions |T5[ are quasi-scaling type. Scaling is 
violated in pionisation region. But this almost not influences properties of families, since 
they are formed mostly by particles from fragmentation region. Scaling is violated a little in 
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fragmentation region as long as the interaction of a primary particle of PGR occurs not with 
a nucleus of air. All these models are in the agreement with accelerator's experiments and 
represent extrapolation of their properties to superhigh energies. They differ insignificantly 
in details of the interaction: in quantity of use flavors of secondary particles, in a way of 
inclusion of diffraction processes and production of jets with large transverse momentum. 
In the majority of quasi-scaling models the intensity of families at normal chemical com- 



position is 2-4 times more then observed [1^, |20|, ^ . A dilemma is arisen: either admit 
scaling violation in the fragmentation region or assume heavy composition of PGR 
This explains the appearance of simulations based on compositions with dominance of heavy 
elements (Table |). By this point of view paper is especially interesting as different pos- 
sible decisions of the dilemma are proposed in it. 

MGO model is used in the present work for the analysis of data observed with the 
help of X-ray emulsion chambers. The Model is based on the theory of quark-glion strings. 
Diffraction processes, generation of jets with large transverse momentum, production of 
strange and charm particles are included in it. The important peculiarity of MGO is a large 
inelasticity coefficient. Details of the model are described in |^. 

Let us note, that for to day the most popular is code GORSIKA with some variants 
of strong interaction model. Gomparisons of MGO predictions with results of GORSIKA 
calculations in variants DPM and VENUS were made in work p2| , . There was shown 

that MGO predicts faster absorption of hadron component than both variants of GORSIKA. 
As it will be apparent from the following fast absorption of hadron component provides the 
agreement of MGO model with observed intensity of 7-families. Model MGO and GORSIKA 
differs also in the average characteristics of families. The characteristics in MGO are closer 
to the experimental data than in GORSIKA. This explains a choice of MGO as base for the 
subsequent analysis. 

With the help of MGO model nuclear-electromagnetic cascades in the Atmosphere gen- 
erated by various nuclei were simulated at power like energy spectrum. Index of integral 
energy spectra of PGR, 7, for all nuclei was also simulated. 7=1.7 for energies less then 
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3 X 10^^ eV and 7=2.2 below the "knee" were input in this case. 

There are 6 sections in the present paper. Besides Introduction (sec.l) and conclusion 
(sec. 6) it contains 4 main parts each of which are subdivided on two subsections concerning 
experimental and modeling considerations. Section 2 is dealing with the problem concerning 
intensity of famihes at the altitude of Pamir. Observed (subsec.2.1) and predicted (sub- 
sec. 2. 2) intensities are considered there. In the next two subsections (3.1 and 3.2) properties 
of experimental and simulated families are described. In subsection 4.1 the consistency of 
the model with the experimental data is discussed while in subsection 4.2 this subject is 
analyzed with the account of possible systematic errors in the experiment. One- and multi- 
dimensional methods of recognition of families generated by iron are elaborated and applied 
to the Pamir data. This is done in subsections 5.1 and 5.2. 5.1 is a short review of the 
situation but 5.2 is an original investigation. In the last section the results are summarized. 

2. Intensity of 7-hadron families . 

In the next two subsections intensity of 7-hadron families are studied from experimental 
and theoretical point of view. In the end of subsection 2.2 the results are compared and 
a solution of the dilemma, heavy composition or scaling violation in the fragmentation 
region, is found. In the frame of quasi-scaling model good agreement between observed and 
predicted intensity at the normal chemical composition is achieved. This is due to large 
inelasticity coefficient in MCO model. 

2.1 Intensity of 7-hadron families, Pamir experiment. 

A group of 7-quanta (true 7 and also e^, e~) and hadrons {ir'^, n", n, p...) produced in 
NEC by interaction of a PGR particle with a nucleus of air is called the 7-hadron family. In 
this paper we analyze families satisfying the following conditions: 

lOOTeV < T.E^ < lOOOTey, 
> 10, E^, El > iTeV, 1cm < < 15cm (2) 
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In (2) E2, is visible energy of a hadron or a 7-quantum, i?^ - average radius of the 
family count off from its energy weighed center. 

Only 7-quanta seated at distance less than 15cm from the energy-weighted centre of a 
family are included in it. 

So-called, aggregation procedure was used at processing of families. It consists of taking 
a pair of 7-quanta (7j7j), being at distance Rij smaller then 0.15 mm, as one with = 
E^. + E^.. The majority of such quanta actually combine their spots of darkness and the 
scanner's eye apprehends them as one. The rest we unite by the aggregation procedure as 
in real so in simulated families. Thus all 7-quanta with R^j <0.15 mm are united. 

Let us pay attention to the condition <lcm. This selection condition is applied to 
families for the first time. It is connected with two circumstances. On the one hand the 
characteristics of families with R^ <lcm are strongly distorted by process of formation of 
dark spots on X-ray film of emulsion chamber. A large part of them are overlapped, some 
of them are almost completely united. These peculiarities of the registration are difficult to 
reproduce in simulated families. Therefore we prefer to exclude them from consideration. 
The same reasons concern to families with TiE^ > 1000 TeV. In their central part there are 
very narrow bunches of 7-quanta which are also difficult to separate. On the other hand 
it was noticed that the relative part of narrow families in experiment is more than 
in simulations (Fig. la) and their properties differ from properties of the other families. 
For example, families with R^ <lcm in experiment are almost away of hadrons and their 
spectrum of E^ is very similar to the spectrum of E^ in purely electromagnetic cascades. It 
also indicates that group of families with R^ <lcm is preferable to exclude from the analysis. 
Figure lb shows that after exclusion of events with R^ <lcm experimental and simulated 
distributions of R-y become closer. 



Intensity of 7-hardon families at Pamir level (4370m above sea, 596 gr/cm?) without 
condition R-^ >lcm is equal to 

(0.69±0.15)m-Var"^str-^ (3) 
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Error bars include statistical and possible systematically uncertainties. 
Total number of families studied in the present work including narrow events {R^ <1 cm) 
equals to 226. Without them it becomes 174. 

Thus the experimental intensity of families satisfying the criterion (2) is 

lexp = (0.69 ± 0.15) X 174/226 = (0.53 ± Q.12)m-^year-^str-^ (4) 

2.2 Intensity of 7-hadron families , MCO-model. 

The vertical intensity of 7-hadron families can be expressed as: 

I^^ = E[Ia{E > Eo) X eA{E > E^) x ^0/^^^] (5) 

Ia{E>Eo) - is an intensity of nuclei with atomic number A and energy E > Eq, 
^a{E > Eq) - is an efficiency of families production satisfying the given conditions. 
In turn 

SA = Nf,jNA{E > Eo) (6) 

Na{E>Eq) is the number of NEC, induced by a primary particle with atomic number A 
and N^°'"^ - number of families, which they have produced. 

It is convenient to have Ia in percentage of I, i.e. to normalize it on total intensity: 

Ca = Ia/I (7) 
In (5) flo and flfam are solid angles for primary particles and families: 

no = 271 nfam = 27l/{l + T/Xatt) (8) 

Here T is air pressure at the altitude of installation's exposition (600 gr/cm^ for Pamir) and 
Xatt is attenuation mean free path of families. 

As one can see from (5-8) at determination of calculated intensity model defines only sa 
(efficiency of production), whereas Ca and the index of energy spectra of nuclei are set by 
assumed chemical composition. Simulation calculations have shown that Xatt is independent 
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of A in the limit of their errors and equals to 75 ± 6 gr cm ^. This value is in a good 



agreement with the experimental figure A^^f = 78 ± 4 grcm ^ JTO 



Vertical intensity for all particles of PGR 

r{E > Eo) = T.I\{E > Eo) (9) 

is estimated up to very high energy. It is customary to use empirical expression for TjI\{E > 
Eq) given in article [Q. We also use it but change the integral energy spectra index -1.6 
accepted in [Q to more recent value 7=-1.7. As a result we have: 
r{E > lO^^eV) = (50 ± 20) m-^year-hrr^ 

The new figure for differential vertical intensity is obtained in [|1|. It gives for PIE > 
10^5): 

r{E > lO^^eV) = (47± 12) m-^year-^srt-\ 

Combining these two results we use the following value for P: 

I{E > lO^^eV) = (50 ± 15) m-^year-^srr^ (10) 

Efficiencies obtained by means of MCO model for different nuclei are brought in Table |T|. 



N^"^{E > Eq) are given in the first row of Table|T|. N^'^{E > Eq) are actual number 
of cascades simulated isotropically over zenithal angle interval from 0° up to 43°. For light 
nuclei (P, He) Eq = 10^^ eV. In the case of groups of nuclei CNO and SiMg and also for 
Fe nuclei the minimal primary energy was put as 

= v4°-^ X lO^^eV (11) 

Such a choice of Qmax and Eq is determined by the fact that nuclei with energy less than 
specified Eq and O > Qmax form no more than 2% from total number of families. 

In the second row of TablePTT] corrected number of cascades Na{E > 10^^ eV) are brought. 

Na{E > lO^^eV) = N'/^/il - cosA3")/{E^y-'^ (12) 
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TABLE III: Efficiency of family production by different nuclei, e^, and average energy of 
nuclei, responsible for them, E^. 





P 


He 


CNO 


SiMg 


Fe 


P(7=-2.2) 




63.6 


87 


30 


20 


30 


62 


Na X 10^ 


237 


324 


912 


1113 


3450 


231 


jsjfam 


682 


416 


461 


434 


756 


409 


EA X 100 


.29 ± .01 


.13± .01 


.051 ± .003 


.038 ± .003 


.022 ± .001 


.18± .01 




15. 


29. 


53. 


66. 


80. 


9. 



In Table |T| N^"""^ is number of families generated by the given nucleus and - efficiency 
of their production. Data for protons having above the "knee" {E^nee = 3 x 10^^ eV) spectral 
index 7 = -2.2 are given in the last column of Table [IT| . 



Figures given in Tables | and |T| allow to determine efficiency of family production, e, at 
the given chemical composition, Ca-, and energy spectrum index, 7. Further on one can get 
predicted vertical intensity of families, JJ^^ , average energy of primary nuclei responsible 
for families, Efam, and average In A. 



e = S(Ca X eA)/m 



A 



(13) 



= P.^Ca X X (1 + T/Ai)] (14) 
Efam = ^{Ca X X EA)/^iCA X e^) (15) 



(InA) = X Inv4)/SCA (16) 

The values Ca and (In A) for investigated CC are given in Table | (section 1), eA and Ea - 
in Table |Tl[ Parameters calculated by (13, 14, 15) are brought in Table Dependencies 
of family production efficiency on In A for pure nuclei and on {In A) for various chemical 
compositions are shown in Figure 2a. Predicted intensities of families at different chemical 
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compositions are brought in Figure 2b. In Figure 2b the strip hmited by broken hues 
corresponds to region allowed by the experiment. Table |I^ and Figure 2b show that only 
normal compositions (NC) with 7 = -1.7 for all nuclei and NC with 7 = - 1.7 for all 
components except protons (7P = - 2.2 for energy more then Eq = 3x 10^^ eV) are consistent 
with the experimental intensity of 7-hadron families. The heavy compositions 0, ^ predict 
intensity of families essentially less then the experimental value. Let us note, that an account 
of possible systematic errors in energy measurement of 7-quanta and hadrons only strengthen 
reliability of the conclusions. Discussions in subsection 4.2 show that predicted intensity 
at normal composition become (0.58 ± 0.18) m^'^year^^ srt~^ and at heavy composition - 
(0.25 ± 0.08) m^'^year~^ srt~^ after the account for systematically errors. 
Thus we come to the following conclusions: 

• Model MCO eliminates the dilemma: heavy chemical composition or strong scaling vi- 
olation in fragmentation region. In framework of quasi-scaling models the agreement 
between experimental and calculated intensity of families at normal chemical compo- 
sition is managed. This progress is due to rather large inelasticity coefficient in MCO 
model. 

• Heavy and superheavy chemical compositions in MCO model give too low values of 
family's intensity. 

One more result of simulations should be underline. Table ^ indicates that the average 
energy of nuclei responsible to family's production is about 10^^ eV. Therefore the conclusion 
following from our investigations of families correspond to the energy interval just above the 
"knee" of energy spectra of PCR. 

3. Characteristics of 7-hadron families. 

7-hadron families are characterized by a number of measured parameters. Their defini- 
tions and descriptions are given in the next two subsections. Again in the first one exper- 
imental questions and problems connected with measurement of parameters are discussed. 
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TABLE IV: Calculated parameters (see the text) 



Composition 


e X 100 


If am 


Efam,PeY 


Normal lH 


0.16 


0.71 ± .22 


22. 


Heavy!, § 


0.076 


0.35 ± .11 


32. 


Superheavy m 


0.051 


0.23 ± .07 


45. 


Normal(7p = -2.2) 


0.11 


0.49 ± .15 


17. 



The results of calculations and some specific circumstances are analyzed in the second sub- 
section. On the base of simulated events parameters sensitive to atomic number of incident 
nucleus are found out. Short conclusions about agreement between experimental and cal- 
culated values of parameters at the normal composition are given in the end of subsection 
3.2. 

3.1 Characteristics of 7-hadron families, Pamir experiment. 

7-hadron families are characterized by number of measuring parameters. Conditionally 
they can be divided into 4 classes: 

• 1. Characteristics of 7-quanta related to energy: n^, HE^, HE^/n^. 

- number of 7-quanta, Si?^ - total energy of 7-quanta. 

• 2. Spatial characteristics of 7-quanta: R^, E^R^, = T^E^R^/ E^ and d. 

R-y - radius of a family (average distance from the center of the family), E^R^ - average 
product EyR^ , R^ average radius weighted by energy. Parameter d = riin/nobs is 
defined as the ratio of the number of initial 7-quanta, nj„, to the number of observed 
7-quanta, Hobs = n^- An initial 7-quantum is responsible for a narrow group of spots 
on an X-ray film, which are the result of an electromagnetic cascade induced by it 
in the Atmosphere. Observed dark spots being on a small distance, Rij, from each 
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TABLE V: Average values of parameters, P, their statistical errors, cxp, in experimental 
families and sensitivity of parameters, S. 





nh 


iiy(cm) 


r>E 




d 


Qn 




p 


3.1 


2.8 


2.4 


27. 


.63 


.11 


.14 


ap 


0.3 


0.1 


0.1 


2.0 


.01 


.01 


.01 


S 


1.38 


1.33 


1.26 


1.04 


.96 


.65 


.60 



other are combined into one initial 7-quantum, if Rij/ {1/Ei + l/Ej) < 10 TeV x mm. 
Using this algorithm, called decascading procedure, the number of initial 7-quanta is 
determined. 

• 3. Characteristics of hadrons related to energy: nh, '^E'l, 
qe = T,Ell{T,E\ +T,El), = nj(n^ + rih) 

Uh - number of hadrons, E^ - energy transferred by a hadron into soft component 
(visible energy of a hadron), E£^^ - total visible energy of hadrons, qe - fraction of 
energy carried by hadron component of a family and qn - fraction of hadrons in total 
multiplicity. 

• 4. Spatial characteristics of hadrons: Rh, E^Rh and etc. Let us note, that the last 

characteristics are not examined in this work, since the number of hadrons 

rule, small and consequently their spatial characteristics have very wide fluctuations. 

Parameters belonging to a given class are subject to common systematic errors. In the 
1-st class they are determined by errors in energy measurement of 7-quanta. For energy from 
4 up to 50 TeV the relative error is about 20%. The effects of saturation of darkness and the 
overlapping of spots appeared for larger energies. Whenever possible these effects are taking 
into account during primary processes of families. At the next stage of treatment, 7-quanta 
being on distance Rij < 0.15 mm are united in one and only families with T,E^ < 1000 TeV 
and R^ > 1cm are included into analysis (see subsection 2.1). As it seems these restrictions 
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eliminate the main part of systematic errors of the 1-st class parameters. This question is 
discussed in more detail in subsection 4.2. Complexities in determination of spatial char- 
acteristics (the 2-nd class) were already discussed in subsection2.1. We should remind that 
aggregation of quanta with Rij < 0.15 mm and exclusion of families with Ry < 1 cm from 
analysis set aside most difficulties. Main uncertainties of the 3-rd class's parameters are 
connected with determination of visible energy of a hadron, E^^, based on its darkness. Due 
to two reason the situation here is easier than with 7-quanta: the saturation of a darkness 
comes later, as distribution of density is more flat in case of hadron and overlapping is not 
present, since average distance between hadrons is much more than between 7-quanta. How- 
ever, the relation between visible energy and density of darkness for hadron are investigated 



not so carefully as for 7-quanta [^. This complexity will be discussed in subsection3.2 in 
connection with determination of visible energy of hadrons in simulated families. Average 
values of parameters of the experimental families with their statistical errors are given in 
Table 0. 

Sensitivity of a parameter to atomic number of primary particle is defined as: 

S = {<Pf,> - <Pp >)/Dp (17) 

< Ppe > - is an average value of the given parameter for families induced by Fe, < Pp > 
- the same for families induced by protons. Dp - dissipation of a parameter for primary 
proton. S was calculated by means of simulated families for primary protons and iron. Let 
us note, that all parameters are defined in such a way that < Pp > is less than < Ppe >■ For 
this purpose in two cases it was necessary to depart from initial definitions of parameters. 
Parameter d was introduced for the first time in work |^ as the ratio d = riots /nm- We 



have redefined it by replacing d — > 1/d. Parameter qe, was introduced in work as 
qe = T,E^/ {T.E'^ , +T,E2). We transformed it as g^; — > 1 — qe- In Table parameters are 
brought in order of decrease of their sensitivity. Characteristics of families with S < 0.5 
are not given there. At a research of chemical composition they can be only harmful. Not 
having sensitivity they are useless but the systematic errors in them can enter distortions 
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into final results. Distributions of two parameters sensitive (R-y) and not sensitive (n^) to 
atomic number of a primary nuclei are demonstrated in Figures 3a and 3b. 



The average value of a given parameter of families, P, at certain chemical composition 
can be expressed by a formula: 



where Pa is an average value of the given parameter in families generated by a nucleus with 
atomic number A; Ca - a fraction of nuclei A at the given chemical composition of PGR; 
- an efficiency of family production by nuclei A with the energy above lO^^eV. As above, the 
model determines values Pa and e^, whereas CC and the power index of energy spectra of 
nuclei A are set as an input of simulations. It was already noted that all procedures used at 
selection and processing of experimental data were applied to simulated families. Peculiarity 
of simulated events is a modelling of hadron registration in X-ray emulsion chamber and 
determination of energy transferred by it into the soft component, E]^ = x Eh. A hadron 
is registered by X-ray film if it has interacted in the chamber and energy transferred by it, 
E!^, is more than threshold value - 4 TeV. The probability of hadron interaction and factor 

is determined by a design of X-ray emulsion chamber. For Pamir carbon chambers 
special investigations have shown that the probability of interaction is about 0.7 and 

has a distribution, f{K^), similar to incomplete 7-function 



For example, if a = 1.5 and (3 = 0.075 average < > has quite reasonable values equal to 
0.188. However, the condition E^ > 4 TeV makes < > dependent on Eh (Figure 4a). 
Therefore the more important parameter is not (Kj) but K^ff 



3.2 Characteristics of 7-hadron families, MCO-model. 




(18) 




(19) 




(20) 
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TABLE VI: Average values of parameters, Pa, and their dissipation, Dp, in simulated 
families. 





nh 


(cm) 


r>E 




d 


qn 


Qe 


p 


2.7 


2.8 


2.4 


24. 


.59 


.10 


.10 


Dp 


2.4 


1.5 


1.6 


15. 


.16 


.07 


.11 


p a 


2.5 


2.7 


2.3 


24. 


.60 


.10 


.10 


Dp 


2.4 


1.4 


1.5 


14. 


.16 


.08 


.10 


He 


3.3 


3.2 


2.9 


30. 


.63 


.12 


.20 




2.8 


1.6 


1.7 


17. 


.16 


.08 


.11 


CNO 


4.3 


4.0 


3.6 


34. 


.69 


.14 


.15 


DcNO 


3.2 


1.7 


1.9 


17. 


.13 


.08 


.10 


SiMg 


5.0 


4.5 


4.1 


38. 


.72 


.16 


.16 


DsiMg 


4.0 


1.8 


2.0 


19. 


.13 


.08 


.11 


Fe 


6.3 


4.8 


4.4 


40. 


.74 


.18 


.19 




4.7 


1.9 


2.1 


18. 


.13 


.09 


.11 



Here F{Eh)dEh is the energy spectrum of hadrons in families. As calculations shown for 
simulated families Keff = 0.23. This value is quite compatible with that offered in p9| . 



Distributions of obtained in |29| for Pamir chamber and our approximation of f{K^) 
are given in Figure 4b. They are in a reasonable agreement. Knowing probability of hadron 
interaction, its energy and distribution of and using Monte-Carlo method one gets for 
each hadron answers to questions: whether hadron interact in the chamber? what energy 
it transfers into soft component = {KjEh)'^ and whether it is registered, i.e. E^ is 
more than Ethr = 4 TeV? Such algorithm was applied to simulated families to get hadron 
characteristics: Uh, q-n, Qe, etc. Average values of sensitive (5* > 0.5) parameters of families, 
Pa, and their dispersions. Dp, for various primary nuclei(P" - Protons with energy spectrum 



having " knee" ) are brought in Table ^ 
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TABLE VII: Expected values of families characteristics at various chemical compositions. 





nh 




Rj 




d 


Qn 


QE 


Normal 


3.2 


3.0 


2.6 


26. 


.62 


.11 


.11 


Heavy 


3.6 


3.3 


2.9 


29. 


.64 


.12 


.12 


Superheavy 


4.2 


3.7 


3.3 


32. 


.66 


.13 


.14 



Dependencies of Uh, R-y and don A are shown in Figures 5a and 5b. As above characteris- 
tics of families for all nuclei are calculated at integral energy spectrum index 7 = —1.7. Data 
for families generated by protons having power spectrum with bend in point i^^o = 3 x 10^^ eV 
are brought in the third and fourth lines. Up to the bend 7 = —1.7, after 7 = —2.2. Table|V^ 
shows that the "knee" of the spectrum insignificant influences on average characteristics of 
proton induced families. 

One can calculate expected values of parameters for the given CC (Table [VII] ) with the 
help of expression (18) and using Tables |, |T| and ^ where various compositions, Ca-, 
efficiency of family's generation, e^, and the average values of parameters. Pa, are given. 

Comparison of Tables |V| and |V11| shows that characteristics of families at normal compo- 
sition are in good agreement with the experimental data whereas predictions for heavy and 
the more so for superheavy compositions differ much from the observations. This subject is 
discussed in the next two subsections. 



4. Comparison of the experimental data with the results of MCO model. 



Detailed comparison of the experimental data with the results of MCO model is given in 
the following two subsections. In the first one test is applied to the above mentioned 
sensitive parameters but before parameters uncorrelated between each others were found out. 
All of them show a much better agreement with the results corresponding to the normal 
chemical composition than that at the heavy compositions. In the next subsection not only 
sensitive but other main parameters of families are attracted to the analyze with the aim 
of study the possible systematic errors in the experiment. Here also good agreement with 
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calculations at the normal chemical composition is achieved. 

4.1 Compatibility of the experimental data with the results of MCO model. 

Before making the final conclusions about chemical composition of PGR by comparison 
of experimental data (Table [V|) with predictions of MCO model (Table |V11| ) one has to find 
out which of the studied characteristics are not correlated. Otherwise traditional (for ex- 
ample, x^) ciiid not traditional (neural net) approaches can incorrectly estimate a degree of 
agreement of experimental and simulated data. Among selected seven parameters sensitive 
to chemical composition (Table |V| ) there are 2 groups of strongly correlated characteristics. 
They are characteristics related to the energy of hadrons, Uh, qn, Qe (for example see Figure 
6a) and spatial characteristics of 7-quanta, R^, , E^R^ (Figure 6b). Stands to apart pa- 
rameter d, which weekly correlate with both of the groups (Figure 6c). Parameters belonging 
to the different classes do not correlate also (Figure 6d). We have chosen three sensitive and 
not correlated parameters: n^, R^, and d. Only they are used in the subsequent analysis. 
For comparisons of the experimental and simulated families the following quantities were 
calculated: 

Xp = [{Pexp — Pmod) / Cf Pexpf (2l) 

and sum for three not correlated parameters 

X3 = ^{nhexp — nhmod) / O'Uhexpf + 
\^P'~iexp R-ymod) / ^^R^exp ]^ + [{dexp - dmod)/o-dexp]'^] / 3 (22) 

Here Pgxp and Pmod are an average value of some parameter but crPexp is an error of Pexp- 
Instead of total error ((Tg^,p + cTmod) only corresponding a^^^ stands in (21) and (22) since 
statistical errors in calculations are much less than experimental uncertainties. The results 



are shown in Table VIII 



Each should be near one if the experiment and calculations are in a good agreement, 
since the number of degrees of freedom for separate parameter is 1 and Xs is an appropriate 
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TABLE VIII: Values Xp for various parameters and xl- 



cc 


Xnh 






xi 


Normal 


0.11 


3.4 


0.95 


1.5 


Heavy 


2.7 


20. 


0.24 


7.7 


Superheavy 


13. 


54. 


7.2 


25. 



with the 3 degrees of freedom divided on 3. This expectation is fulfilled only for the 
normal composition, but is not satisfied for the heavy and superheavy compositions. It is 
necessary to notice, that rather large values of Xit-y =3.4 for the normal composition can 
indicate on the presence of some systematic errors. This problem will be analyzed in the 
next subsection. The dependence of x| on < InA > is shown in Figure 7. Dotted line 
corresponds to xl =4. Confidence level at =4 is less than 1%. Figure 7 shows that all 
compositions with < InA > more than 2.5 are above the line and therefore can be excluded 
with confidence level more than 99%. 



4.2 Consents of MCO model vi^ith the experimental data. Systematic errors. 

It was shown in the previous section that MCO model with normal chemical composition of 
PCR is in a reasonable agreement with experimental data. Nevertheless two circumstances 
cause doubt. Firstly, experimental average radius of families = 2.8 ±0.1 appears to 
be less than calculated radiuses at any chemical compositions = 3.0 - 3.7 (Table |V| 
and |VII| ). This is the cause of a little bit large values Xi?,7=3.4 and accordingly X3=l-5 
for the normal CC (Table |V1I1| ). For heavy compositions these quantities are much more. 
Secondly, experimental average values of = 21 ± 2 is less than calculated = 24 — 26 
at any chemical composition. Till now was not interesting as it does not depend on 
atomic number of incident nuclei (see Figure 3b). There is a suspicion that R^ and 
are subject of some systematic experimental errors and it is desirable to find their common 
origin. In this subsection we want to carry out the more complete analysis of the consent of 
MCO model with experimental data considering possible systematic errors at measurement 
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TABLE IX: Values of parameters, P, and their Xp- 



Exp and CC 


Par am. 


rih 




d 


q-n 






E'y 


Exp 


P 


3.1 


2.8 


.63 


.11 


27. 


21. 


12.1 




ap 


0.3 


.01 


.01 


.01 


2. 


2. 


0.5 


Normal 


P 


3.2 


3.0 


.62 


.11 


26. 


24. 


11.1 




xl 


0.11 


3.4 


.95 


.06 


1.3 


13. 


3.1 


Heavy 


P 


3.6 


3.3 


.64 


.12 


29. 


25. 


10.8 




xl 


2.71 


20. 


.24 


1.5 


0.86 


20. 


7.0 


Superheavy 


p 


4.2 


3.7 


.66 


.13 


32. 


26. 


10.5 




xl 


13. 


54. 


7.2 


9.5 


5.2 


25. 


12.4 



of parameters. Values of the most important parameters characterizing 7-hadron families 
measured in the our experiment and predicted by MCO model for various CC together with 
xl are brought in Table |X[ As already was noted the most serious deviation of model 
from experimental data is displayed in parameters and n^. For heavy composition this 
disagreement is even stronger than at the normal CC. The common reason of all large 
deviations can be systematic errors in energy determination, E^, by measurement of spot's 
darkness, D, on a X-ray film. The algorithm of transition from D to is rather complex. 
Besides theoretical connection between energy of 7-quantum and the number of electrons 
in a circle of certain radius under given thickness of lead where a film is located (transition 
from E^ to ng) it includes the property of the film (ng —>■ D curve, saturation of D at large 
rie etc), condition of the exposition of the film (its background, gap between the film and 
lead plate) and so on. All these effects mostly are accounted and controlled. However, it is 
not excluded that some of discrepancy can nevertheless be left. Group effects in families are 
added to these problems - partial overlapping of spots and even their aggregation (the last 
effect we took into account at processing of simulated families). To study an influence of 
energy measurement errors on family's characteristics we have added to simulation of events 
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modeling of distortion by replacing Etme to E, 



meas ■ 



For this we used some distortion functions, f{E): 




'true 



(23) 



Two distortion functions were tested : 



fl{Etrue) = W/Etrue X (Etrue/^O) 



a 



(24) 



with a >1 (a <1) up to 10 TeV and a <1 (a >1) after 10 TeV to provide underestimation 
(overestimation) of true energy 



with /? <1 in the whole interval of Etrue- Parameters a and (3 were varied. Typical exam- 
ples of distortion functions are shown in Figure 8. At the particular chose of a, function 
fi{E) reflects underestimation of energy up to and after Etrue=^^ TeV with different depen- 
dence of Emeas oii Etrue- Oppositc to fi{E), f2{E) Overestimates true energy up to 70 TeV 
and then underestimates it. f2{E) was taken from work[^. An attempt of modeling of 
passage of 7-quanta through the lead plates of the chamber and then its fixing on a film 
with the accounts of peculiarities mentioned above was carried out there. Various values 
of a were set in calculations to ensure 10%, 20%, 50%, (-20% and -40%) underestimation 
(overestimation) of energy near measured threshold 4 TeV and —5%, —10%, —20% under- 
estimation of energy at true energy lOOTeV. Dependencies of Xp for the normal chemical 
composition on a degree of energy underestimation or overestimation in Figure 9a ,9b. The 
Figures show that all Xp have quite admissible values already at 10% of underestimation 
of energy {Etrue = 4.4 TeV instead Emeas= 4 TeV) and no systematic errors at lOOTeV 
{Emeas = Etrue at Etrue >100 TcV). It docsu't meaus abscucc of systematic errors at large 
energy Etrue >100 TeV. Simply it appears that values Xp weekly depend on an underes- 
timation of energy at Etrue >100 TeV (Figure 9b). As the 10% systematic error is quite 



and 



f2{Etrue) = 70/ Etrue X {Etrue/70) 



(25) 
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TABLE X: Xp values for normal and heavy chemical compositions at 10% and 20% under- 
estimation of energy near the threshold. 



cc 


Etrue 






d 


qn 








xi 


Normal 


4.4 


0.09 


1.85 


.11 


.02 


0.11 


3.4 


1.37 


0.70 


Heavy 


4.8 


0.63 


7.39 


2.65 


1.84 


1.88 


4.63 


2.15 


3.60 



possible we conclude once more that MCO model at the normal chemical composition is 
in a very good agreement with the experimental data. Situation at heavy and the more 
so at superheavy compositions is different. In the previous section it was shown that they 
contradict to the experiment. Figure 9c, similar to 9a, for heavy composition shows that 
no one distortion function can improve the consent of MCO model with the experiment at 
heavy composition. Xp for normal and heavy chemical compositions at the most favorable 
distortion functions for each of them are given in Table We should note that in the case of 
superheavy composition Xp values are so large that this composition is out of consideration. 

Finally let us discuss the predicted intensity of families at used distortion functions. 
Since after introduction of distortion functions the number of events, satisfying the selection 
criteria of families, decreases, diminished also expected intensity of families. For the normal 
composition at Etme = 4.4 TeV the intensity become equal to Ifam = (0.58 ± 0.18) m~^ 
year~^ srt~^ . For the heavy composition at Etme = 4.8 TeV Ifam = (0.25 ± 0.08) m~^ 
year~^ srt^^. Let us remind that experimental intensity is (0.43 ± 0.12) m^^ year~^ srt~^. 
From the given figures it is apparent that after corrections for possible systematic errors 
the conclusions made in section 2 concerning chemical compositions based on intensity of 
families strengthen. 

We have used a number of procedures simulating process of 7-quanta and hadrons reg- 
istration in X-ray chamber. These are: a procedure of aggregation (unification near-by 
spots); an ascription to each hadron a factor of energy transfer into soft component, Kj, 
distributed according to incomplete 7-function; possible underestimation of energy threshold 
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of 7-quanta and hadrons, Ethr', an account for systematic underestimation of energy of the 
most energetic 7-quanta (replacement in simulations Efrue to Emeas)- This can appear to 
be not correct and artificial. Actually, as any complex installation, X-ray emulsion chamber 
brings in certain distortions in NEC falling on it. As a result instead of true parameters, 
Ptrue, one has to do with their measured values, Pmeas- To have an opportunity to compare 
parameters measured in experiment with that of simulated, it is necessary to model pro- 
cesses occurring in X-ray emulsion chamber. In Pamir collaboration the routine, carrying 



out this problem is named the code of "passages through the chamber" pO[]. This code is 



rather complex and is an analogue of routine JEANT PTI , used at accelerators for imitation 



of distortions, connected with this or that installation: ATLAS, CMS and others p2|, p3 



In the present work "the program of passage through the chamber" was divided on stapes 
(aggregation, i^'^, Etme E^eas) to look after influence of each of the effects separately. By 
such approach we have resulted in the conclusions about a role of each type of distortions 
in determination of measured parameters of families. 

5. Selections of families generated by iron nuclei. 

The attempt of selections of families similar to that induced iron(Fe)nuclei by means 
of image recognition methods were done in the several papers. All of them claim that 
compositions enriched with heavy components are in contradiction to experimental results. 
Only compositions near to the normal can satisfy them. The description of the method and 
short review of the mentioned works are given in the subsection 5.1. In the next subsection 
our original investigations are brought. The most encouraging is the result of reconstruction 
the fraction of Fe-like families inputted into calculations by the method, which is proposed. 
The main conclusion is again only results at the normal composition are consistent with 
the experiment. 

5.1 Selections of families generated by iron nuclei. A short review. 

The attempt to determine a fraction of iron nuclei in PCR on a database of X-ray-emulsion 



chamber was undertaken in work [27, 34, 35] for the first time. Almost simultaneously by 
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two groups of authors it was offered to apply to family multidimensional analysis of an image 
recognition. The sense of it is that using a difference in distributions of some parameter of 
families generated by iron and proton (see Fig. 3a) one defines a limiting value for the given 
parameter, Pum- Event with Pfam > Pum, is considered as a family similar to one generated 
by iron. Such an event is attributed to group of families "Fe". Otherwise if Pjam < Piim, 
then family is added to group "P". A fraction of families, /', " similar to iron " is thus 
defined. Multidimensional method of an image recognition means that the limiting values 
are defined for several parameters P^m, Pfimi Pum ■■■■ Families which simultaneously satisfy 
the requirements: 

pi > pi p2 > p2 

are attributed to group " Fe" . 

Fraction of families induced by nuclei A and satisfying the limiting conditions ( Fe-like 
produced events) are designated Ra- 

Two types of errors in the image recognition method. The error of the first type is to 
attribute families from proton to "Fe" group {Rp ) and the error of the second type is not 
to "recognise" iron family and not to added it to "Fe" group, (1 — Rpe)- The quality of 
selection is defined by values of these two errors. As more sensitive are the parameters to 
atomic number, the stronger differ distributions of parameters for P and Fe families. As a 
consequence cleaner and more complete is the selection as both errors are small. The error 
of the first type includes all families from other primary nuclei (except Fe) falsely attributed 
to the "Fe" group. 

Quasi-scaling models were used for determination of limiting values of parameters Pum * 
in three mentioned works. Training sets of families generated by P and Fe were simulated 
and distributions, similar to that shown in Figure 3, were obtained. These distributions 



allowed determination of the boundary quantities Pum *. It was found Q that it is 
impossible to use simultaneously more than two parameters because of limited statistics of 
families. Therefore either E^R^ and asymmetry parameter of family, b, or E^R^ and were 
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used in [Q. Another parameter of asymmetry, a, and parameter used in the previous 



sections, or their combination with were involved in |^. Parameters d, a and in 
different combinations were analyzed in [^. 

The fraction of families attributed to group "Fe", /', is equal to 

f' = fxRF, + {l-f)x Rp (26) 

where / is the real fraction of families generated by iron. From (26) 

f = U'-Rp)/{RFe + Rp) (27) 



Authors of the all three works ([^, 0, |3^) applying the described method to families of 



Pamir Collaboration have found that the fraction of families generated by iron does not 
exceed 2-3 %, i.e. that the compositions of PGR with domination of iron contradict to the 
experiment. 

The method of selection families generated by iron was essentially improved in work [§, 
where experimental data of Pamir-Chacaltaya Collaborations were used. In contrast to us 
(see conditions (2)) in papers based on Joint Pamir-Chacaltaya experiment a group of 7- 
quanta and hadrons with TjE^ + S-E^ >100 TeV at and El > A TeV was considered 
as a family 0. For recognition of an image of a family produced by iron author of [6!!] 
used a neural net method, with the help of which mult i- dimensional analysis is reduced to 
one-dimensional. 15 parameters describing a family were used in 0: n^, n^, Si?^, T^E^ and 
so on. With the help of neural net the set of them was reduced to one, ?/p, and condition 
which attributes a family to the "Fe" group was yp >0.5. The quality of selection making 
by this way has appeared to be rather high: Rp ^ (1 — Rpe) ~ 15%. 

Also the next step on the way of CC study was made in [^]: not only families produced 
by P and Fe but also generated by other nuclei (He and CNO, SiMg groups) were examined 
there. It means that these nuclei also contribute to group "Fe". Following our consideration 
in this case(24) should be transformed to: 

/' = X Ra) (28) 
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and accordingly 

/=(/'- S'(/^ X RA))/RFe (29) 

Here is a fraction of families generated by nuclei A, Ra - their fraction faulty attributed 
to group "Fe" (except Rpei Rpe is true portion), /' - the part of families satisfying the 
limiting condition. In S(/a x Ra) contributions of all nuclei are summarized, in S'(/a x Ra) 
iron families are not included. Let us underline once more that we designate true value of 
Fe induced families as fpe and estimated by (28) as /. The training sets define Ra and the 
given chemical composition of PGR determines fA (see (1) and Table f|. 

The conclusion of is the same as in previous works ^ : the heavy composition 



of PGR is excluded by experimental data concerning 7-hadron families. 

For the sake of completeness of the review refer to work In this work rather complex 
parameters, reflecting structural properties of families, and intensity of " structural" families 
are used. The conclusion in is opposite to the above-stated. Authors affirm that the 
superheavy composition of PGR is necessary for description of considered properties. As 
one of the critical remarks to ||^ let us note that authors of ||^ did not show that the 
model they used describes the more simple characteristics of families at the superheavy 
composition. 

5.2 Selections of families generated by iron nuclei. Original consideration. 



We have used methods similar to mult i- dimensional and one-dimensional 
approaches as a following step of chemical composition researches. Before starting GG 
analysis we have found parameters n^, i?^, and d sensitive to A and are not correlated 
between each other (section 3). In contrast to |]^, ^ we used only them. The parameters 
are analyzed in their reduced form: 



Xp = {P-P,)/5Pp 



(30) 
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Here P is a value of some parameter in a given family either experimental or simulated, 
Pp - average values of the same parameter in families induced by protons, SPp - disipation 
of this parameter in proton families. 

In reduced form of variables all distributions are dimensionless and for families generated 
by protons are dispersed around average values Xp = 0. At multidimensional analysis 
simultaneously three conditions were used to attribute a family to "Fe" group: 

^rih > limi > Xji^ lim-i > X^ Hm (31) 

For one- dimensional analysis a new parameter X3, was introduced: 

X3 = {Xn,+XR^+Xd)/3 (32) 

For this parameter a value X3 nm was also found and the family was attributed to group 
"Fe" if its 

> ^3 lim (33) 

Limiting conditions Pum were determined with the help of integral distributions of Xp for 
families induced by P and Fe using training sets of simulations. Integral distributions at once 
show what fraction of proton families (Rp) at the given Pum will be by mistake attributed 
to "Fe" group and what fraction of iron families will not recognized as Fe-like, (l-Rpe)- An 
example of such distribution is given in Figure 10. The arrows in Figure 10 indicate limiting 
values of X3. One of them corresponds to about equal errors -Rp = (1 — Rpe) ~ 20%, while 
the other - to Rp = 5% and (1 - Rpe) ~ 40%. 

It is apparent that both errors Rp and 1 — Rpe depend on chosen limits um, X^^ n^, 
Xd lim, and X3 lim- We investigated two sets of limiting parameters. First, at which Rp = 
5%, i.e. only 5% of proton families are falsely attributed to "Fe" group, and second, such 
that the errors of the first and second types are approximately equal -Rp (1 — -Rpe)- True 
values of fraction iron induced families, /pe(see Table IH), and expected fraction of families. 



selected to group "Fe", /', for different chemical compositions are brought in Table XI at 
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TABLE XI: True fraction of Fe families, fpe, and expected fraction of Fe-like families,/', at 
various chemical compositions 





CC 


fFe%, true 


/'%, one-dim. 


/'%, multi-dim. 




Normal 


2.8 


13. 


11. 


Rp = 5% 


Heavy 


16.4 


21. 


19. 




Superheavy 


31. 


30. 


28. 




Exp. 




12. ±3. 


11. ±3. 




Normal 


2.8 


26. 


27. 


Rp = 20% 


Heavy 


16.4 


35. 


37. 




Superheavy 


31. 


47. 


47. 




Exp. 




29. ±4. 


25.±4. 



two sets of limiting parameters Xp Let us remind that fp^ is fully determined by the 
given CC while /' is the result of processing of simulated data for set of limiting parameters. 

Table ^ is composed using training sets for various nuclei P, He, CNO, SiMg and Fe. Ra 
was determined for each nucleus and then the sum /' = x R^) was found. We should 
like to pay attention to the fact that for normal composition the true iron families figural 
spiking is lost among families selected in "Fe" group. At the normal composition the true 
Fe fraction fp^ = 2.8% while Fe-like proton /' varies from 11%-to 27% depending on Rp. 

To check efficiency of the used method to determine a part of Fe induced families we 
have investigated how it reproduces fpe at various chemical compositions. The results (/) 



are given in Table PCII| for one particular selection (multi-dimensional selection, Rp = 20%). 
Results of the other three selections ( mult i- dimensional Rp = 5%, one-dimensional at Rp = 
20% and 5%) are identical. 

In Table [XII rows correspond to simulated compositions, columns to compositions by 
means of which corrections were done using (29) . For comparison / with the true values of 
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TABLE XII: Fractions of iron induced families, / , in chemical composition of families, 
calculated by (29). Mult i- dimensional selection, Rp = 20% 





P 


Normal 


Heavy 


Superheavy 


Fe 


M,% 


p 


0.2±2.5 


-4.3± 


-2.6± 


-1.2± 


27.±2.5 





Normal 


8.7±1.7 


4.2±1.7 


5.9±1.7 


7.3±1.7 


36.±1.7 


2.8 


Heavy 


21. ±2. 4 


16. ±2.4 


18. ±2.4 


19.±2.4 


48.± 


16.4 


Superheavy 


35. ±2.7 


30.±2.7 


32. ±2. 7 


34.±2.7 


62. ±2. 7 


31. 


Fe 


79.±6.4 


75.±6.4 


76.±6.4 


78.±6.4 


106. ±6.4 


100 


Exp. 


5.9±5.0 


1.4±5.0 


3.1±5.0 


4.5±5.0 


33.±5. 


? 



fpe the letter for each simulated compositions are given in the last column of the Table. 



Table pCII| demonstrates rather encouraging result. Independently to the composition 
used for corrections (for exception of pure Fe), the fractions / are close to corresponding 
true value. As it is expected the best agreement / with true fpe is obtained if the composition 
for corrections is close to the "real". The corresponding figures are underlined in Table |X11 . 
From the above the following procedure of processing of experimental data is suggested. 
After determining an experimental value /' corrections should be done for various CC, for 
example for five compositions testing in Table Xll . Then receiving preliminary result (five 
values for /) one takes that at which chemical composition determining corrections is closest 



to obtained /. This is shown in the last row of Table pCII| . The final value of found fraction 
of Fe induced families is underlined. 

The total analysis of our experimental data was as follows. Four values of fractions of 
families similar to iron, /', were found corresponding to two type of the methods (multi- 
dimensional and one-dimensional) and two sets of limiting parameters for Rp = 5% and 
Rp = 20%. They are brought in the last row of Table ^l[ Values / were reproduced 
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by the help (29). Quantities / are given in the last row of Table PCI1| corresponding to 
various correction compositions. Table PCI1| shows, that only normal composition gives self 



consistent values /. In this case true values fpe = 2.8% and set of experimental quantities 
of / fluctuates from 1.4% to 4.5%. In the case of the heavy composition the experimental 
fraction / = 3.1% contradicts to true value fpe = 16.4%. Even more disagreement shows 
the superheavy composition: true fpe = 31% while obtained f is equal to 4,5%. 

As it was specified above we have tested four selections criteria of Fe-like families. Ta- 
ble pCll| is given for one of them. The results of all other selections are identical: events. 



generated by iron nuclei, compose about (2 - 3)%. This is in agreement with the normal 
composition and sharply contradicts to compositions enriched by iron. 

6. Conclusion 

The purpose of the present work was an investigation of Chemical Composition of Primary 
Cosmic Rays with energy around 10^^ eV close to the "knee" of energy spectrum of PCR. In 
this region of energy information about CC of PCR can be obtained only by an indirect way 
on the base of study of extensive air showers or families of 7-quanta and hadrons registered 
by X-ray emulsion chambers. 

In the present work the data of Pamir and Pamir-Chacaltaya Collaborations concerning 
families were used. The analysis of a material was made by means of comparison of the 
experimental families with that simulated by quasi-scaling model MCO 

First of all it was shown that MCO at normal CC (close to composition at lO^^eV, about 
40% P and 20% Fe) predicts intensity of families in complete agreement with experimental 
observations (section 2). Thus long-term dispute, heavy chemical composition or strong 
scaling violation in the fragmentation region was solved. 

Further it was shown that not only intensity, but also all main characteristic of families 
(they are about 15) are well described by MCO at normal chemical composition (section 3). 

Characteristics of families sensitive to atomic number of an incident nucleus were found 
out on the base of simulated events generated by primary protons and iron. Among them 3 
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parameters not correlating with each others were chosen. They are number of hadrons Uh, 
radius of a family and parameter d, describing electromagnetic structure of an event. 
It has appeared that their average values in experiment are close to the same in artificial 
families at the normal composition of PGR and are in disagreement with the heavy (15% 
P and 58% Fe) and, especially, with the superheavy (7% P and 70% Fe) compositions, xi 
was calculated for the three above mentioned parameters. In the case of the normal CC it 
is equal to 1.9, for two iron enriched CC X3= 7.7 and 25. correspondingly (subsection4.1). 

Despite of the rather good agreement of experiment and calculations at the normal com- 
position two basic characteristics of families and have rather significant deviations 
from expected values: Xn^ = 13. and X/?^= 3-4- The investigations have shown that if one ad- 
mits 10% underestimation of energy near to registration threshold of 7-quanta (about 4 TeV) 
and introduces the appropriate distortion function into simulated families, then correspond- 
ing become equal to Xn^= 3-4 and x|;^ = l-85 (subsection4.2). At 20% underestimation 
of energy for all characteristics of families are close to 1. On the other hand any under- 
estimation or overestimation of energy do not bring to an agreement between experimental 
average characteristics of families to those simulated at heavy and superheavy compositions. 

Another attempt to investigate chemical composition was based on the selection of fam- 
ilies similar to that generated by iron. Such method was already applied to experimental 
data of Pamir [^, 31, ^ and Pamir-Chacaltaya Collaborations 0.1n all four investigations 
fraction of Fe-generation 7-hadron families was small and equal 2-3% . Thus it appear that 
the heavy compositions contradict to the fraction of families from Fe in the experimental 
data. 

In this article we continued elaboration of an application of image recognition methods to 
CC research by means of families (sections 5). Limiting values {Xp um) for the three above 
mentioned sensitive parameters (P) in their reduced form Xp were determined (section 5). 
Families with Xp > Xp were attributed to "Fe" group. A united parameter X3 = 
{Xn^ + Xfi^ + Xd)/S with its X3 lijn was also used. Ability of the method was investigated 
by the help of simulated families. Sets of families generated by pure protons and irons as 
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well as by the normal, heavy and superheavy chemical compositions were analyzed. Families 
"similar" to iron were selected for all these sets. The fractions of Fe-like families and the 
restored ((29)) parts of "true" Fe families were estimated. For each set of simulations the 
found "true" fraction of Fe-families were in coincidence with the fraction of Fe families in 
the given chemical composition. 

The method of selection of Fe-like families was applied to experimental events. A frac- 
tion of families similar to that induced by iron was found. It appeared to be close to the 
appropriate fraction at the normal composition. After corrections by means of (29) the part 
of families from iron in the experimental set constitutes (2 - 3)%. Such value is in complete 
agreement with the normal CC and in sharp contradiction to the iron enriched compositions. 

Thus the basic conclusions of the given work are: 

MCO model at the normal chemical composition completely agrees with the experimental 
data of 7-hadron families; 

The chemical composition of Primary Cosmic Ray in the energy region near to 10^^ eV 
just above the "knee" of its energy spectrum is close to the chemical composition at energy 
around 10^^ eV; 

The models of nuclear electromagnetic cascade in the Atmosphere with chemical compo- 
sitions enriched by heavy elements contradict to the experimental data on families. They 
predict too low intensity of families, incorrect values of the characteristics of families and 
too large fraction of families generated by iron. 

One more conclusion very important for the particle physics follows from the present work. 
All the quasi-scahng models [|, 0, |T^, |T5|, ^ |2^, except MCO, give 2-3 times larger 
intensity of families than the observed one. It means that in these models hadron absorption 
mean free part is too large to provide proper intensity of families. Smaller mean free part 
in MCO is due to large an inelasticity coefficient than in the other models. The better 

is 



agreement of average characteristics of families in MCO than in other models ||22|, ^ 
also the result of stronger absorption of hadrons. Therefore large an inelasticity coefficient 
is necessary feature of inelastic interactions at superhigh energies. 
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FIG. 1: Rj - distributions. Continues line - experiment, dashed - normal CC. a) all families, 
b) families with >lcm. 




FIG. 2: a) Families production efficiency as function of In A. o - pure nuclei, A - different 
CC. b) Intensity of 7-famihes at different CC. 
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FIG. 3: Distributions of (a) and (b) for families generated by protons (continues line) 
and iron (dashed line). 




FIG. 4: a) Dependence of (K-y) on E^. b) K^ff distribution. Histogram is taken from p9| , 
curve is used in this work. 
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FIG. 6: Correlation: a)bctween rth and q^^ k=0.60±0.02; b)between and 
ERy, k-:0.88±0.01; c) between d and R^, k-0.34±0.03; d) between rih and R^, 
k=0.00±0.04. 
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FIG. 7: Dependence of xl on (In A). 




FIG. 8: Two types of distortion functions. 
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FIG. 9: Dependencies of Xp and Xs on: a) degree of energy underestimation near measured 
threshold 4 TeV at normal chemical composition, b) degree of energy overestimation near 
true energy 100 TeV at normal chemical composition , c) the same as a) for heavy CC. 
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FIG. 10: Integral distributions of for families, generated by protons and iron. Limited 
values of correspondent to conditions Rp — Rpe and Rp — 5%, are shown by dashed 
lines. 



